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ABSTRACT 

Context. The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. 
Aims. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ 
Helioseismic and Magnetic Imager (SDO/HMI) data. 

Methods. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version 
of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 
arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured 
from HMI white light data. 

Results. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the 
two components of the horizontal velocity v x and v y measured in the sky plane and the simultaneous line of sight component from 
SDO/HMI dopplergrams Vp, we derive the spherical velocity components (v r , v„, v 9 ). The azimuthal component v 9 gives the solar 
differential rotation with a high precision (+0.037 km s~') from a temporal sequence of only three hours. 

Conclusions. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial 
and temporal resolutions from hours to months and years with the SDO data. 

t***~ ' Key words. The Sun: Atmosphere - The Sun: Granulation - The Sun: Convection 

CN 

CO 1. Introduction of the Sun in white light at low (SDO/HMI) and high (Hinode) 

' spatial resolution. 

CO ' The dynamics of the Sun's surface is one of the major elements 

t-H | in understanding the time evolution of its magnetic activity. It is 
^ • a real challenge to measure the surface motions at all spatial and 
temporal scales and compare them to those coming from the sim- 
ulations. Recently, the HMI instrument aboard the SDO satellite 
allowed us a new step in that direction. By following the proper 
motions of the solar granules, representative of solar plasma evo- 
lution, it is possible to define the flow field on the solar surface 
dRoudier et al.ll2012h from a small spatial scale of 2.5 Mm up 

to nearly 85% of the solar radius (Fig. 6 of that paper). We mea- In this paper, we describe a comparison between velocity 
sured the velocities in a cartesian coordinate system, where x and fields projected on the sky plane that are obtained in the south 
y denote the coordinates in the sky plane with x parallel to the pole r egion with Hinode data and SDO/HMI da ta using the CST 
direction of the solar rotation and z directed towards the observer code dRieutord et al. 2007t iRoudier et al.ll2012l) . In the next see- 
along the line of sight. Beyond 0.85 R , the v x and v y compo- tion we discuss the corrections done in order to get an accurate 
nents appeared to be noisy, but the v x component showed a trend velocity close to the solar limb. From the v x , v y and Doppler ob- 
indicative of the solar rotation. In order to identify whether im- servations, we describe the transformation to the local surface 
provements could be made in the determination of the horizontal velocities v r , v g , v<p. Finally, we present an application to mea- 
velocities beyond that limit, we used simultaneous observations sure solar differential rotation with a short time sequence (3h) 
up to high latitudes with low noise. Discussion and conclusions 
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Fig. 1. v x computed from LCT and CST on Hinode observations Fig. 2. Angular rotation (sidereal) computed from LCT and CST 
on December 10, 201 1. on Hinode observation on December 10, 201 1. 



2. Observations 

2.1. Hinode observations 

We used data sets from the Solar Optic al Telescope ( SOT), 
onboard the Hinod e Q mission (e.g., ISuematsu et alJ [2008; 
Ichimotoetalj[2004). The SOT has a 50 cm primary mirror with 
a spatial resolution of about 0.2"at a wavelength of 550 nm. For 
our study, we used blue continuum observations at 450.45 nm. 
from the Hinode/SOT BFI (Broadband Filter Imager). The ob- 
servations were recorded continuously on December 10, 2011, 
from 16:11:33 to 19:07:05 UT. To get the limb as reference, 
the south pole was observed at the position reported in the 
Flexible Image Transport System (FITS) header of Hinode, i.e, 
Xcem = -10.08" and T cent = -969.84". The time step was 45 
sec and the field of view was 1 1 1 .6 " x 1 1 1 .6 " with a pixel size 
of "1089. After alignment, the useful field of view reduced to 
104" x 82". To remove the effects of the oscillations, we applied 
a subsonic Fourier filter. This filter was defined by a cone in k-co 
space, where k and u> are spatial and temporal frequencies. All 
Fourier components such that oj/k > V C ut-off = 7kms~ werere - 
moved so as to keep only convective motions dTitle et alj l989). 

2.2. SDO/HMI observations 



The HMI dScherrer et ai]|20Tl ISchou et alj|2012b onboard the 
Solar SDO provides uninterrupted observations over the entire 
disk. This gives a unique opportunity for mapping surface flows 
on various scales (spatial and temporal). Using the SDO/HMI 
white light data on December 10, 2011, from 16:11:15 to 
19:06:45 UT, we derived horizontal velocity fields from im- 
age granulatio n tracking using the newly developed version of 
the CST code (Ro udier et alj|2012l) . The time step was 45 sec- 
onds with a pixel size of 0'.'5. The solar differential rotation dis- 
cussed in Sect. 6 was determined from SDO/HMI white light and 
Doppler data taken on August 30, 2010 from 8:00:45 to 1 1 :09:45 
UT. 



1 Launched in 2006 the Hinode spacecraft , was designed and is now 
operated by JAXA (Japan Aerospace Exploration Agency) in coopera- 
tion with NASA (National Aeronautics and Space Administration) and 
ESA (European Space Agency). 



3. Hinode south pole flow field 

Measuring the solar rotation near the poles is a difficult task be- 
cause there are few solar structures at high latitudes suitable for 
tracking. However, the Hinode images allow us to follow so- 
lar granules with good contrast up to the solar limb. First, we 
compare results of two methods to measure the solar motion 
close to the south pole using Hinode observations: local cor- 
relation tracking (LCT) ([November & Simonl [l988b and CST 
dRieutord et alJl2007t iRoudier et al.ll2012l) . Both methods track 
the horizontal motions of granules in the field of view. More pre- 
cisely, the LCT method of obtaining the horizontal velocity uses 
a spatial window that simultaneously accounts for the solar gran- 
ules and integranular structures and may cover several granules. 
In contrast, the CST method measures the velocities by follow- 
ing the trajectory of each granule, i.e., solar plasma, during the 
life of the coherent object, which is defined by its appearance and 
disappearance if the granule does not split or merge. The data 
were aligned and the P angle evolution corrected by 0.018°/hour 
in order to get a perfect co-alignment with SDO data used in the 
following. On December 10, 2011, the Bq angle was -0.25°, 
resulting in insignificant projection effects with a negligible evo- 
lution during the observation period. Figure[TJshows good agree- 
ment of the v x component from the LCT and CST up to latitude 
75°. Close to the limb, the effect of the spatial window used in 
the LCT prevents correct velocity determination. Figure Ogives 
a solar sidereal rotation rate close to the pole of around 1 1 °/day, 
which corresponds to a polar period of 32.73 days . This is in 
agreement with previous determinations (Beck 2000). 



4. Comparison of the velocities from SDO and 
Hinode 

Our main goal is to compare the surface flows measured with 
a high spatial resolution (Hinode, 0.1") and a low resolution 
(SDO, 0.5"). We extracted the same field of view from the SDO 
data set as the Hinode one and performed a very precise co- 
alignment of both sequences. To get the best co-alignment, the 
south polar region was observed and the solar limb was taken as 
an absolute reference. An additional check was performed after 
the alignment process by locating the brightest features (facular 
points) at the beginning and end of both sequences (Hinode and 
SDO). A very good match of the structures indicated very good 
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Fig. 4. v x component from Hinode and SDO Hinode (0.1' 
SDO (0.5"), same field and duration. 
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Fig. 5. Vy component from Hinode (0.1") and SDO (0.5") same 
field and duration. 



alignment of both sequences during the 3 hours. We then ap- 
plied the CST to both sequences to get the horizontal velocities 
(v x and v y ) over the entire field of view. Figures |4] and |6]show 
a good agreement up to 78°of southern latitude. This indicates 
that SDO observations can be used to determine the solar rota- 
tion in the central meridian region up to significantly high lati- 
tudes. However, as shown in Fig. [5] the meridional component 
measured on SDO data exhibits an offset of 0.4 km s relative 
to the Hinode component, which was used as the reference (be- 
cause the velocity is close to zero at the limb). Figure [5] shows a 
good correspondence between Hinode and SDO velocities when 
the 0.4 km s _1 offset is removed (curve marked as SDO-0.4). 
We tried to elucidate the origin of that offset by dete rmining all 
possible errors as listed in Table II of IStrousI £000). The esti- 
mation on the SDO measurement gives a total error on the v y 
component of around 0.034 km s , which is ten times smaller 
than the measured offset. 

As described above, one of the major differences between 
the two sequences is the pixel size, which is 0.1" and 0.5" for 
Hinode and SDO, respectively. In order to analyze the sequences 
under the same conditions, we degraded the pixel size of the 
Hinode observation to the SDO one: 0.5". We then applied CST 
to the Hinode degraded sequence to get the horizontal veloc- 



ities, which were also k-a> filtered. Figure [7] still shows good 
agreement for the v x component and, in the same way, for the 
siderial angular rotation, at the latitudes up to 78°. The offset of 
0.4 km s observed in Vy component of the Hinode data at 0.5" 
is clearly visible in Fig. [7] We conclude that the offset is caused 
by the combination of the lower spatial resolution and decreas- 
ing contrast to the limb. One has to bear in mind that we only 
observe the south pole region in detail, where v y is practically 
identical to the radial component of the flow. Thus, our observa- 
tion can be generalized so that the radial component is always 
affected by the offset and we can correct for it. Indeed, different 
processes play a role in generating that radial component: at very 
high heliocentric angles, one pixel covers several granules. The 
granule part close to the limb is of lower contrast: it tends to be 
lost in the segmented images. This introduces an artificial radial 
motion of some granules that we see as the offset in the v y com- 
ponent on the central meridian when comparing high-resolution 
(Hinode) and low-resolution (SDO) measurements. This radial 
effect was previously observed in our measurements, but its ori- 
gin was not identified. Since we now have the origin of that error, 
we can measure it and correct for it all over the Sun. One way to 
get the radial offset is to average the radial velocity component 
over a circle centered on the solar disk. Due to the Bo angle, the 
radial correction is not necessarily symmetrical in the northern 
and southern hemispheres of the Sun. This is why we treat the 
average process in the northern and southern regions of the Sun 
separately. Figure|8]shows, for example, the plot of the measured 
average radial component in the southern part and the overplot- 
ted fit obtained by using a polynomial function of the fifth de- 
gree. Figure [9] shows the entire profile of the correction of the 
radial velocity to be applied to the velocities v x and v y all over 
the Sun (north and south), allowing a good representation of the 
velocity field around 90% of the solar radius. In the following, 
this correction is applied systematically to all velocity fields. 



5. Determination of the spherical components of 
the velocity 

From the velocity components v x and v y measured in the sky 
plane and the simultaneous line of sight velocity \'o measured 
from the SDO/HMI dopplergrams (e.g., Fig. [TTJ, we can de- 
rive the spherical velocity components v r , vg, v v , projected onto 
spherical coordinates r, 9, <p using 



v r (9, <p) = cos 9 * sin tp * v x 

+ (sin 9 * cos Bq - cos 9 * cos tp * sin Bq) * v y 

+ (cos 9 * cos if * cos Bq + sin 9 * sin Bq) * vd op 

vq{9, (p) = - sin 9 * sin ip * v x 

+ (sin 9 * cos (p * sin Bq + cos 9 * cos Bq) * v y 

+ (cos 9 * sin Bq - sin 9 * cos <p * cos Bq) * vd op 

v 9 (9, ip) = cos ip* v x 

+ sin tp * sin Bq * v y 

- sin tp * cos Bo * vdo P 

The system of coordinates is also shown in Fig. [10] We note 
that tp is along the longitude and 9 is along the latitude. The co- 
alignment of all components has been checked carefully by us- 
ing the information from FITS headers of the SDO data. Before 
using the dopplergrams, some data reduction had to be applied. 
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Fig.3. SDO (left) and Hinode (right) observations on December 10, 2011, at 16:11:15 and 16:11:33 respectively. The location of 
the Hinode field is shown on the SDO image 




First, the dopplergrams were averaged over the same time inter- 
val as the v x and v y sequence, her e three hours. The limbshift had 
to be be corrected (Ulrich 2010) and, in this case, we used the 
limbshift function determined by the SDO team over the month 
of August 2010 (P. Scherrer, private communication). 

By defining z — 1 - cos p, where p is heliocentric angle, 
the limb-shift correction (in km s _I ) is given by limb-shift(z) = 
-0.664z + 0.775z 2 + 0.284z 3 . The mean velocity of the central 
region of the dopplergrams was taken as zero origin. The dopp- 
lergrams were resampled to the same size as the v x and v y com- 
ponents. 

To convert the v x and v y velocities in the sidereal system v^ ld 
and v^. ld , the Earth's orbital displacement was taken into account. 
This correction depends on the P angle as 

v^ ld = v x + v e cos P 



Vy = v y + v e sin P 

On August 30, 2010, it was 0.13 km s" 1 on v x and 0.04 km s" 1 
on v y , where v e is the Earth velocity component projected on the 
Sun's surface. 

Figure [12] shows the resulting velocities v r , vg, v v for the 
three-hour sequence on August 30, 2010. The v r map exhibits the 
radial component where the downflow is visible in the sunspot 
regions. Close to disk center we observe a lower contrast of v r , 
due to the low contrast in Doppler velocity in this part of the 
Sun because of the mostly horizontal flow on the Sun's surface. 
The v, component is not well determined close to the limb be- 
cause the projection effects are predominant. The component 
clearly shows the solar differential rotation with a lower velocity 
close to the pole. To our knowledge, this is the first time that we 
can measure and visualize the solar differential rotation with a 
three-hour time sequence. The latitudinal component vg map al- 



4 



Roudier et al.: Comparison of solar horizontal velocity maps measured using SDO/HMI and Hinode data 







Fig. 11. 



voop velocity maps for the three-hour sequence on August 30,2010, and the Intensity at the beginning of the sequence. 



lows us to essentially observe motion on supergranulation scales 
with a comparatively short time sequence. 



6. Determination of the solar differential rotation 

One of the first scientific ap plications of the CS T algorithm on 
SDO/HMI data described in lRoudier et al.1 d2012h was to deter- 
mine the solar rotation from the granule displacements. From the 
longitudinal velocity component v v , it is possible to calculate the 
solar differential rotation. In order to reduce the noise, we aver- 
age v v over bands in longitudes. Figure [13] shows the computed 
differential rotation averaged over bands limited to longitude of 
(-2,2), (-10,10), (-20,20), (-40,40) degrees. The overplot- 
ted continuous line repres ents the solar rotation me asured by 
the spectroscopic method dHoward & Harvevlll970l) . which is 
our reference to evaluate the noise of our measurements. The 
noise level is found to be 0.94, 0.37, 0.27, 0.267 °/day for the 
set of above-mentioned longitudinal bands. Thus, the average 



equatorial solar rotation determined for these four bands around 
the central meridian is 1.99 km s" 1 ± (0.133, 0.052, 0.038, 
0.037) kms~' for the different longitudinal bands. 

Our results show that, for the first time, using the CST we can 
get a determination of the solar rotation with a time sequence as 
short as three hours at very high precision, namely, 1.9% in the 
best case. 



7. Discussion and conclusion 

The comparison between high and low spatial resolution obser- 
vations from the Hinode and SDO satellites allowed us to quan- 
tify the quality of the horizontal velocities determined over the 
full Sun. We found that the differential rotation along the central 
meridian is well determined up to latitude 60°. With a longer 
time series, higher latitudes will be reached on depending the 
Bo angle. Comparison of the meridional components shows an 
offset of 0.4 km s _1 due to a combination of three factors: the 
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Fig. 13. Solar differential rotation on August 30, 2010, computed for the four different angles: 2°(top left), 10°(top right), 20°(bottom 
left), 40°(bottom right) either side of the central meridian. 



low spatial resolution, the limb gradient contrast, and the seg- 
mentation process. However, we describe a way to measure and 
correct for the radial effect all over the Sun and get nearly a full 
Sun velocity measurement. From the velocities v x and v y mea- 
sured in the sky plane and the simultaneous line of sight velocity 
from SDO/HMI dopplergrams, we derived the spherical veloc- 
ity components (v r , vg, v^,). From the longitudinal component, it 
is possible to get the solar differential rotation with high preci- 
sion ( +0.037 km s _1 ) using a temporal sequence of only three 
hours. This is remarkable because other methods require at least 
one month of data. That particularity opens a new field of study 
of the solar rotation and motions over the solar surface. In this 
way, we can revisit the dynamics of the solar surface at high spa- 
tial and temporal resolution from hours to months and years with 
the SDO data. In particular, it will be of great interest to com- 
pare CST convective velocities wi th numerical simula tions and 
constraints from helioseismology dGizon & Birchll20Tl . 
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Fig. 6. Angular solar rotation from Hinode and SDO Hinode 
(0.1") and SDO (0.5"), same field and duration. 
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Fig. 8. Determination of the radial offset close to the south solar 
limb and the adjusted function 
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Fig. 9. Entire profile of the correction of the radial velocity to be 
applied on the velocities v x and v y all over the Sun 
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Fig. 7. The comparison of v x , angular solar rotation and v y and 
from Hinode (0.1"), Hinode (0.5") and SDO (0.5") 



7 



Roudier et al.: Comparison of solar horizontal velocity maps measured using SDO/HMI and Hinode data 




